Collective dynamics of interacting Ising spins: Exact results for the Bethe lattice 
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We study the low temperature dynamics in films made of molecular magnets, i. e. crystals 
composed of molecules having large electronic spin S in their ground state. The electronic spin 
dynamics is mediated by coupling to a nuclear spin bath; this coupling allows transitions for a 
small fraction of electronic spins between their two energy minima, 5*^ = ±S, under resonant 
conditions when the change of the Zeeman energy in magnetic dipolar field of other electronic spins 
is compensated by interaction with nuclear spins. Transitions of resonant spins can result in opening 
or closing resonances in their neighbors leading to the collective dynamics at sufficiently large density 
Po of resonant spins. We formulate and solve the equivalent dynamic percolation problem for the 
Bethe lattice (BL) of spins interacting with z neighbors and find that depending on the density of 
resonant spins Po and the number of neighbors z the system has either one (2 < 2 < 6) or two 
{z > 6) kinetic transitions at Po — Pci ~ e^^^'"^ /{3z) and Po = Pc2 ~ e^^/z. The former transition 
is continuous and associated with the formation of an infinite cluster of coupled resonant spins 
similarly to the static percolation transition occurring at Po ~ The latter transition, z > 5, is 
discontinuous and associated with the instantaneous increase in the density of resonant spins from 
the small value ~ 1/z to near unity. Experimental implications of our results are discussed. 

PACS numbers: 7080.Le, 72.20.Ee, 72.25.-b, 87.14.Gg 



I. INTRODUCTION 

Percolation theory describes the flow in heterogeneous 
media. It is successfully applied to a variety of physi- 
cal, chemical, biological and even social processes ranging 
from the hopping conductivity in doped semiconductorsi 
to the evolution of large genetic networks and self- 
organized criticality2i^. For instance, the low tempera- 
ture conductivity of hopping insulators can be modeled 
by the equivalent network of resistances replacing ele- 
mentary electron hopping eventsi. The percolation the- 
ory finds the optimum set of interconnected open chan- 
nels characterized by resistances not exceeding a certain 
maximum and forming an infinite cluster to guide elec- 
trons through the sample. 

Open percolation channels are usually treated as sta- 
tionary in time. It is relatively easy to study the static 
case numerically. Moreover there exist analytical solu- 
tions including percolation on the Bethe lattice (BL) and 
in low-dimensional systems^'^'i'^. In all cases the perco- 
lation kinetic transition is found to be continuous. Sys- 
tem parameters including, for instance, diffusion coeffi- 
cient, correlation radius and average size of percolating 
clusters all show scaling behavior near the critical point 
similarly to the one for second order phase transitions. 

The model of static percolation is often an approxi- 
mation. Flowing particles interact with each other and 
their motion can open or close percolation channels lead- 
ing to, for instance, conductivity noised. The significance 



of fluctuations for the cooperative dynamics of quantum 
defects, spins in spin glasses, protons in ionic conductors 
and electrons in hopping insulators was pointed out in 
Refs .-'^^i^^'^^i^^'^'^i^^ . In a more refined analysis of perco- 
lation one should include the possibility that the flows 
through open channels can open or close other channels. 
These dynamic correlations can affect kinetic transitions 
in a fundamental way and lead to the cooperative dy- 
namics, e. g. percolation may occur at smaller density of 
open channels and become discontinuous. In this paper 
we suggest a simple model of dynamic percolation for 
the low-temperature kinetics of the system of interact- 
ing magnetic molecules and obtain an exact description 
of the kinetic transitions on BL. According to the argu- 
ments of Refs.— the use of BL is more justified in our 
case than for the static percolation because of the simi- 
larity in the phase space structure for Bethe lattice and 
many-body system of interacting spins. In both cases the 
size of the phase space grows exponentially with the sys- 
tem size (2^ for N spins and z{z — for the number 
of nodes in BL with the coordination number z and N 
shells of the tree). 

In the present work we study the dynamic percolation 
model on the Bethe lattice and then project our results 
on the electronic spin dynamics in films made of magnetic 
molecules such as Mni2, Fes, etc. These molecules are 
usually composed of transition metals, forming the mag- 
netic core, coupled to various organic ligandsi^. They 
possess large electronic spin S in the ground state (both 
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Mni2 and Fes molecules are characterized by the ground 
state "central spin" S = 10). At low temperatures only 
the two lowest states of each molecule, characterized by 
the spin projections S'f = ±5 to the easy axis z, are 
occupied and each molecule can be modeled as a two- 
level system with two states corresponding to the Ising 
pseudo-spin 1 /2. The tunneling gap 2 Aq between the two 
lowest spin states in zero external transverse field is tiny, 
^10~^^ — 10~^Jir. It is much smaller than other relevant 
parameters including the strength of the dipolar interac- 
tion between the nearest-neighbor spins, Ud ~ 0.1 K, 
and the half-width of the distributions of the nuclear bi- 
ases, Eo ~ 10~^ — 10~'^K. Under these conditions only 
spins exposed to a total longitudinal bias (external field 
plus internal dipolar demagnetization field) smaller than 
Eq can efficiently transfer between their two states due to 
energy exchange with the nuclear spin bath^Si^ii^. These 
spins are the resonant, or "open", ones. At the same 
time, since Uo > Eg, each spin flip affects all its neigh- 
bors and may create new resonances (or destroy existing 
ones, see Sec. |IT1). The associated collective dynamics of 
spins is studied within the dynamic percolation model in 
Sees. Uni IV] (see Fig. [I]). 

To be more specific, one may be interested in the cor- 
relation function C{t) = N-^ E,(5f (i)5f (0)) where N is 
the total number of spins, and the average is taken over 
initial conditions, and evolution histories. An interesting 
quantity to investigate in the infinite time limit would 
be the fraction of spin which never change their magne- 
tization Wf. This fraction can be easily seen to result 
in the finite value of C(< — > oo) since dynamically frozen 
spins contribute unity to the sum above. The other ques- 
tion concerns spin diffusion and the possibility to transfer 
energy and magnetization over large distances in which 
case one should be looking for the onset of percolating 
clusters of mobile spins. 

The effect of dynamic and percolation transitions on 
spin relaxation in films of magnetic molecules, based 
on examples of Mni2 and Feg, is discussed in Sec. El 
Our goal here is to investigate the case of completely 
demagnetized samples at relatively high temperatures 
ksT > Ud to avoid the "dipolar ordering" effects and re- 
lated sample-geometry dependent demagnetization fields. 
In this limit the resonant spins are nearly uniformly 
and randomly distributed in the system bulk contrary 
to the case of strongly polarized samples where resonant 
spins form spatially ordered structures, the "resonant 
surfaces" , and initial conditions play an important role 
in the system kinetics^^. 

We predict the reduction of the relaxation rate by 
many orders of magnitude in the vicinity of the tran- 
sition point. At the same time the abrupt change of the 
relaxation rate is smeared out by its continuous depen- 
dence on the total longitudinal bias, which may obscure 
the experimental observation of the transition point. 



II. COLLECTIVE DYNAMICS OF 
INTERACTING SPINS 

Since the model of spin dynamics is already formulated 
in great detail in Refsj ^°i^^'^^ , here we briefly outline its 
main features and then introduce its Bethe Lattice ver- 
sion. As it was mentioned in the Introduction, the weak- 
ness of spin tunneling in molecular magnets allows an 
accurate representation of the system in terms of Ising 
pseudospins 1/2 coupled to each other by the long-range 
magnetic dipolar interaction UijS^Sj with 

1 _ Q„2 p 

t/,, =C/z,a^^^, n,, =^, (1) 

ij 

where C/^i is the strength of the dipole-dipole interaction 
of spins; a is the lattice constant; and Ry is the vector 
connecting molecules i and j. For the sake of simplic- 
ity we limit our consideration to the case of 2D square 
lattice. In this model each spin Si is subjected to the 
"molecular" field of all other spins 

u., = Y,a,,s] (2) 

and the external magnetic field iih, where /i is the related 
magnetic moment. 

The spin dynamics is associated with the spin tun- 
neling between states = ±1/2. The transition of each 
spin between these two states requires them to acquire or 
release the longitudinal energy Eq. ^ to some thermal 
bath because of the energy conservation. The coupling of 
spins to phonons usually responsible for low-temperature 
dynamics in dielectrics is extremely weak and can be ne- 
glected. The relaxation takes place due to the interac- 
tion of the central electronic spins with the nuclear spin 
bath^Si^ii^^. Since in the limit Aq << Eo only spins 
inside the "resonance window" 2Eo are allowed to tun- 
nel, the approximate constraint | fih + Ui |< Eg deter- 
mines the subspace of "resonant spins" . In this form the 
constraint is most suitable for the dynamic percolation 
model studied below in Sections IIIII and IIVI More ac- 
curate dependence of the relaxation time on energy is 
considered in Section |Vl 

The only spins that are allowed to tunnel in any given 
configuration are the resonant ones. Transitions of reso- 
nant spins can change the status of their neighbors from 
non-resonant to resonant and vice versa by changing their 
local bias from being larger than Eo to being smaller 
(see Fig. [TJ At the same time in the new configura- 
tion the transitions of neighbors, diffused into the reso- 
nance window, can push the former resonant spins out of 
resonance. Depending on the density of resonant spins 
Pq ^ Eo/Ud there exist two possibihties. If Pq is very 
small (i.e., the distance between the resonant spins is 
large), the transitions of resonant spins can essentially af- 
fect only their local environment and spins in the sample 
can be separated into a small "mobile" and large "immo- 
bile" groups. The mobile group consists of percolating 



3 



vtJ'ittt^tv 



I A j I t * t t t ^ 



(1) 



(2) 



tween spins. Indeed, such correlations should affect the 
kinetic transition. This assumption is not restrictive 
since, for instance, in Fes and Mni2 the regime of pure 
ground-state tunneling can be reached at temperatures 
T < 0.4 K and T < 0.8 K respectively which is larger 
than Ud in both systemsi^. We also consider only com- 
pletely depolarized system here. 
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FIG. 1: Evolution of "percolating" (or "open") spins affecting 
their neighbors in 2D lattice. When some of the open spins 
(shown by dashed red arrows) make transitions (shadowed 
arrows for steps (2) and (3)) they switch the state of the 
neighboring spins between "open" and "close". 



spins, capable of making transitions in the course of sys- 
tem evolution, while the immobile group consists of non- 
percolating spins which cannot flip in spite of the field 
fluctuations caused by mobile spins. If the density Pq is 
large, the resonant spin transitions create or destroy res- 
onances around them leading to the collective dynamics 
involving practically all spins after sufficiently long time. 

In what follows we investigate the kinetic transition 
between these two regimes. The study of realistic 2D 
lattice of spins, coupled by the dipolar interaction, can 
not be done exactly analytically and requires numerical 
simulations. However, we can solve a similar problem 
on the Bethe lattice with random interactions between 
the neighbors. We assume that spins 1/2 occupy all sites 
of the Bethe lattice and each spin interacts with all its 
z neighbors. The interaction constant Uij between two 
neighboring spins is assumed to be the Gaussian random 
variable with the zero average and dispersion Ud- This 
assumption is approximate - dipolar interactions in a real 
system involve not only nearest neighbors but also spins 
separated from the given one by more than one lattice 
constant. The long range nature of dipolar interactions 
can not be neglected in three-dimensions, but in the two- 
dimensional system the lattice sum of 1/r'^ is finite and 
the nearest neighbor approximation is justified. It is clear 
that the model we study and the real problem differ from 
each other; however they can belong to the same univer- 
sality class because the phase-space structure for inter- 
acting Ising spins is similar to that for the Bethe lattice. 
The phase space grows exponentially with the number N 
of spins as in many-body problemi^iii, which permits us 
to neglect "close loops" for different evolution paths of 
the system in its phase space thus making it similar to 
the Bethe lattice. 

As mentioned above, we restrict our considerations 
to the limit of low temperatures, where each molecule 
can still be modeled as a two-level system, but we as- 
sume ksT > Ud to avoid the dipolar ordering and, 
consequently, any long-range statistical correlations be- 



lli. DYNAMIC PERCOLATION PROBLEM ON 
BETHE LATTICE 

In the Bethe lattice the spin is resonant (open) if its 
local bias (or Zeeman splitting) Ui ~ J2j ^ij^j^ Elq. 
is smaller than Eg- Since the dipolar fields distribution 
in the depolarized limit is similar to the Gaussian one 
(see FigU] below), the fraction of resonant spins (i.e., 
the probability that the spin is open) is determined by 
Pq « 2Eo/(UdV2ttz) < 1. We expect (and this will 
be confirmed by the final result) that the collective dy- 
namics can take place if each spin has approximately one 
resonant neighbor (cf. Refj^), which leads to the esti- 
mate for the transition point as Pqc ^ 1/z. This hap- 
pens at Eo ~ U D \fT^T^\2z) . Near the transition point 
one can approximately ignore correlations between dif- 
ferent configurations of neighbors corresponding to open 
states of the given spin. Indeed, if the spin is open in the 
given configuration the probability that it will be open 
again after k < z turns of neighbors can be estimated as 
Pfe ^ 2Eo/{U D\/2nk). The number of neighbor transi- 
tions Z needed to bring the given spin back to the open 
state can be estimated setting J2k=i ~ 1- One can see 
that at Eo ~ Ud \/7r/(2z) the next opening of the given 
spin can be expected, on average, after turns of around 
z/2 neighbors so the new resonant environment is fully 
different from the previous one. 

Under above assumptions one can formulate the fol- 
lowing model of dynamic percolation on the BL, see Fig. 
[T] Consider a lattice occupied by spins 1/2 having ran- 
dom projection = ±1/2. Assume that in any given 
configuration of spins the rules for deciding about the 
spin states are the same: the probability that a given 
spin is open (i.e., "dynamic") in each of the 2^ configu- 
rations of neighbors is Pq. In other words, the dynamics 
of the model is fixed by having a table for each spin with 
2^ open and close states (the fraction of closed states is 
Wo = I — Pq). Correlations between open states of neigh- 
bors in different configurations are ignored. Only open 
spins are allowed to change their states. As time evolves, 
all open spins have equal chances to make a transition. 
The spin overturn affects all its z neighbors in a way that 
they have to change "status" according to their tables. 
Our main goal is to study the cooperative dynamics of 
spins, characterized by the fraction = 1 — W^,, which 
were involved into dynamics at some stage (if the spin 
ever flips in the course of evolution we call it "percolat- 
ing" ) and to flnd out when an infinite percolating cluster 
is formed by those percolating spins. 
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The problem under consideration has common fea- 
tures and differences with the bootstrap percolation 
problem2^i2£ used as a mathematically idealized model 
for such phenomena as nucleation and growth applied 
in the study of crack formation, magnetic alloys, hy- 
drogen mixtures, and computer storage arrays. It has 
been extensively studied both for the 2D lattice and BL. 
The bootstrap percolation problem can be formulated in 
terms of spins 1/2 placed in all lattice sites. The spin is 
open and allowed to flip if the number of neighbors with 
S^- = +1/2 exceeds the predefined number < n < z. 
The cooperative dynamics is determined by the initial 
density Pq of spins with positive projections and in the 
most investigated case of n = 2 this dynamics vanishes 
in a discontinuous manner at sufRciently small Pq] for 
the 2D lattice the threshold density Pqc approaches zero 
with increasing the system size. In spite of similarities, 
our model differs qualitatively because the bootstrap per- 
colation results in the irreversible evolution of spin con- 
figuration towards the increase in the density of spins 
with positive projection, while in our case the density of 
open spins practically does not fluctuate and we study 
a reversible equilibrium dynamics. One consequence of 
this difference is vanishing Pqc in the bootstrap perco- 
lation problem in the 2D lattice while in our case there 
is no cooperative dynamics at sufficiently small density 
of open spins Pq < 2^^/z (z = 4 in 2D) because non- 
percolating spins form an infinite cluster blocking such 
dynamics, cf. Refii. 

The exact description of the dynamic percolation for 
BL can be obtained similarly to Refsi^i^ (see Fig. [5]). 
First, we calculate the density of percolating spins. 
Note that percolating spins can be in the close state for 
some time. Thus percolating spins include open spins 
and all other spins which can enter an open state at some 
point in time. In what follows we show that the density 
P* can undergo a sharp raise (see Fig. [3]) with increasing 
the density of open spins Pq above some critical value 
Pc2 . This discontinuous transition happens after the for- 
mation of the infinite percolating cluster at Pq = Pd 
which promotes the cooperative dynamics. 



IV. KINETICS TRANSITIONS IN BETHE 
LATTICE 

Consider the probability = 1 — P* that a given spin 
is non-percolating, i.e. it is never involved into dynamics 
despite some of its z neighbors making transitions on BL, 
see Fig. [5] This probability depends on the states of all 
its neighbors, which can be treated as un-correlated due 
to the specific properties of BL. Each neighbor of the se- 
lected non-percolating spin is characterized by the condi- 
tional probability We > W» that it is also not percolating 
(see Fig. [2]). All neighbors are in identical situation by 
construction. 

The probability that the spin is not percolating 
can be determined by considering different local environ- 



ments distinguished by the number of percolating neigh- 
bors A; = 0, 1, ...z, see Fig. [51 There are z\ / {k\(z — k)\) in- 
dependent ways to be surrounded by k percolating spins. 
In each of them the selected spin will experience all 2*^ 
possible states in the course of evolution (assuming that 
percolating spins flip an unlimited number of times; this 
is definitely true before the infinite cluster if formed) and 
in each state it has to remain closed. The corresponding 
probability is (Wq)^ (recall that Wq = 1 — Pq). Summing 
up all the probabilities for different k we get 



E 



z!(l - WefW^ 
k\{z~k)\ 



(3) 



The probability We for neighbors of the given non- 
percolating spin to be non-percolating as well is defined 
similarly to Eq. ^ (see Fig. [2|). One should consider 
only the remaining z — 1 neighbors characterized by the 
same probability We to be non-percolating. Accordingly 
we obtain the self-consistent equation for the probability 
We in the form 



We = 



E 

k=0 



(z-l)!(l-l^e)''I^e 



-k-1 



k\{z-k-l) 



-Wr 



2k 



(4) 



In Fig. Othc solution of Eqs. Q, (g]) for the density 
of percolating spins P* = 1 — W^, is shown for BL with 
z = 4, 5, 6 and 7 neighbors. One can see that for z = 6 
and z = 7 there exists the discontinuous transition at 
Pq = Pc2 (label Pd is reserved for the formation of the 
infinite percolating cluster, see below) where the density 
of percolating spins P* jumps from the small value P* ^ 
1/z almost to unity P* ~ 1 - (Wq)^". At Pq > Pc2 the 
vast majority of spins belong to the infinite percolating 
cluster. 




FIG. 2: Configurations of neighboring spins in the Bethe 
lattice and their probabilities. Red dotted arrows show per- 
colating spins and black solid arrows show non-percolating 
spins. Blue lines connect neighboring spins. 

The physics behind the jump of the percolating spin 
density at the transition point Pc2 can be illustrated in 
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the limit of large number of neighbors z ^ 1 when Eq. 
dH) can be simplified by taking advantage of Pg, Pq ^ 1 
and expanding (1 — Pq)^ w 1 — 2*''Po (recall that Pe = 

1- We) 

Pe « Po(l + PeV-^ « Po CXp(zPe). (5) 

This equation has a solution only for Pq < Pc2 ~ e^^/z, 
which agrees with the numerical solutions for Pc2 at 
z > 7. In our opinion the discontinuous transition is 
caused by an avalanche-type growth of the percolating 
site number in the vicinity of other percolating sites. For 
instance, if the density of percolating sites P* reaches 
the value 1/z then approximately each site has one open 
neighbor. Its probability to become percolating increases 
by the factor of 2 which leads to the formation of new 
open sites. This process leads to the jump in the density 
of percolating spins to near unity. 




FIG. 3: Dependence of the percolating sites density P* on the 
density of open spins Po for z — 4, 5, 6 and 7 neighbors. There 
is a discontinuous change in P* for z > 5 where its density in- 
stantaneously jumps to ~ 1. Transition points Pci ~ 0.0546, 
Pc2 ~ 0.1085 are shown for z = 6. 

The discontinuous transition was found in the boot- 
strap percolation problem on Bethc lattice^^, but at ef- 
fective parameter Pq = l/(z — 1)^. It is interesting that in 
spite of the absence of kinetic transition in the bootstrap 
percolation problem in the infinite 2D lattice, the tran- 
sition in the finite system is essentially discontinuous^. 
This gives us some hope that the discontinuous dynamic 
percolation transition may take place in our problem ap- 
plied to realistic lattices of finite dimensions (rf > 1) sim- 
ilarly to a Bethe lattice. 

The cooperative dynamics in the ensemble of inter- 
acting spins arises in the presence of the infinite clus- 
ter of percolating sites. Such cluster obviously exists at 
Po > Pc2, but it may be formed at smaller density Pq of 
open sites. Indeed, the probability to have a percolating 
site near another percolating site is about a factor of 2 
larger than in the absence of percolating neighbors be- 
cause of the doubling in the number of explored configu- 



rations. The corresponding critical point Pq = Pd where 
the infinite cluster of percolating spins is formed can also 
be found exactly. It can be shown (see Appendix) that 
the formation of the infinite cluster is determined by the 
equation 

(2-l)(l + Po-2Pi) = 
1 + iz- 1)2[P2(1 - Po) + Fi +Fo- 2Pi] , (6) 

where Pm, see Eq. (|13[) in the Appendix, is the proba- 
bility of finding a non-percolating spin surrounded by z 
neighbors where m = 0, 1, 2 of them are definitely perco- 
lating, 2 — TO neighbors arc definitely non-percolating and 
the remaining group of z — 2 neighbors can contain both 
percolating and non-percolating spins. For small Pq we 
have P„ « 1 - 2™Po • (to(1 - We) + W^Y^^. The numer- 
ical analysis of the transition point Pd corresponding to 
the formation of an infinite percolating cluster confirms 
that Pel < Pc2- In the z 3> 1 limit one finds that the 
infinite cluster is formed at Po > Pd ~ e~^/'^/(3z) < Pc2 
in agreement with the numerical solution of Eq. ([5]). 



V. DISCUSSION. HOW CAN THE KINETIC 
TRANSITION BE OBSERVED IN MOLECULAR 
MAGNETS? 

Below we consider the kinetic transition in a 2D square 
lattice of spins representing magnetic molecules using the 
results obtained for the Bethe lattice. This considera- 
tion should be applicable to the recently synthesized two- 
dimensional crystals of Mni2 molecules^. Assume that 
the easy axis is perpendicular to the sample plane. Since 
only spins inside the resonance window can change their 
state, the density of resonant spins in the longitudinal 
magnetic field h can be estimated as 

Po « 2gi^,h)Eo, (7) 

where g{E) is the probability density for the Zeeman 
splitting E (longitudinal bias) originated from the spin- 
spin interaction (see Fig. 3]) and fj. is the magnetic mo- 
ment of the molecule. 

To study the kinetic transition using the previous re- 
sults for the Bethe lattice one has to introduce the num- 
ber of neighbors parameter z. We assume that the num- 
ber of neighboring spins can be estimated as the number 
of spins whose resonance can be affected by the transi- 
tion of the given spin. For the crude estimate we count 
neighboring spins as those coupled to the given spin by 
the interaction exceeding the width of the resonant win- 
dow 2Po- In the limit Eo <€i Ud one can estimate this 
number in a quasi-continuum approach as 

z«7r(C/i,/2Po)2/3. (8) 

According to the solution for the Bethe lattice (see Fig. 
[3]) the sharp change of transition rate for almost all spins 
should take place near the point Pq = Pc2. Assuming 
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FIG. 4: Rescaled probability density for the spin energy split- 
ting on a model square lattice of magnetic molecules. The 
predicted transition point is shown by solid lines. 



z ^ 1 the crude estimate for the transition point can 



be made using the approximate relationship Pc2 
(see Sec. lIVp . Accordingly we get 



Ud 



2{TTeg{tih)UD)^ 



7z 



(9) 



The distribution of the dipolar bias energies g{E) in a 
sample can be easily calculated, see Fig. 21 The mini- 
mum threshold is realized at zero field (maximum den- 
sity of resonances) where g{0) w 0.19/C/c so that we get 

Using the above analysis of the kinetic transition one 
can approximately characterize the spin relaxation for 
different molecular magnets. In Fes one has Ud ~ 130 
mK and i?o ~ 6 mK^^ so even in the absence of the 
external field the system is in the localization regime. 
The situation is different in Mni2, where Ud ~ 70 mK 
and Eo ~ 80 mK— and the collective spin dynamics ex- 
ists at zero magnetic field according to the theory, Eq. 
([9]). However, application of the external longitudinal 
magnetic field can reduce the resonance probability and 
result in the localization. Using Eq. ^ one can esti- 
mate the value of the external field corresponding to the 
transition point from 



g{tihc)UD 



1 

7re 



Up 
2Eo 



1/3 



0.09 



(10) 



and the transition in Mni2 takes place at the field he ~ 
2.8J7d/^ ~ 0.2 T (see Fig|4]). One can also expect that 
the transition rate should decrease near the transition 
point. 

Although theory predicts the absence of cooperative 
dynamics at h > the reality is more complicated be- 
cause spins which do not belong to the transition window 
yet have small but finite transition rate. Indeed, spins 



with the Zeeman energy E > Eg can make transitions; 
however their transition rate is becoming exponentially 



r-i cx exp(- I E \ /£„), 



(11) 



because such transition require simultaneous fiips of large 
number of nuclear spins. Then, we predict the exponen- 
tially small transition rate which can be estimated fol- 
lowing the standard percolation theory approach^ as the 
rate Eq. (fTT|) at the Zeeman splitting E equal to the 
threshold value E^, Eq. ^ 

cx eM-Ec/Eo) « eM-{9{}^h,)/g[^Lh)f). (12) 

This dependence is illustrated in Fig. [5] to show the ex- 
pected strong change in the rate of spin relaxation at 
large fields. Unfortunately, the dynamics slowing down 
takes place at fields h k, 2hc, where the density of res- 
onant spins is already very small (see Fig. [J]) and the 
predicted reduction of the relaxation rate by many or- 
ders of magnitude is more difficult to study. The dy- 
namic transition itself becomes a crossover because the 
transition rate is a continuous non-vanishing function of 
energy, Eq. ([TT|) . 



1x10 



1x10 



1x10 



1x10"' 

2 

FIG. 5: Dependence of the relaxation rate r^^ on the mag- 
netic field near the transition point calculated using the expo- 
nential approach Eq. (|12p . To is the characteristic relaxation 
time. 



In addition to the application of the longitudinal mag- 
netic field one can also affect the transition by diluting 
Fcs magnets. The reduction of the concentration of Feg 
molecules to x ^ 1 will amount for changing Ud/Eq 
ratio, and thus provide a knob for determining the tran- 
sition point experimentally. Assuming that the effective 
constant for magnetic dipolar interaction scales propor- 
tionally to the a;'^/^ one can estimate that at a; < 0.25 
the collective dynamics in Feg films will take place. 
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VI. CONCLUSION 

We considered the dynamic percolation problem on 
the Bethe lattice of spins 1/2. Dynamics was intro- 
duced through the open spins capable to change their 
states and affect the status of neighboring spins switch- 
ing them between open and closed states. The prob- 
lem was solved exactly. We found two kinetic transitions 
including the continuous transition associated with the 
formation of the infinite percolating cluster and the dis- 
continuous transition associated with the avalanche- type 
growth in the number of percolating sites. This model ap- 
proximately describes the low-temperature dynamics of 
molecular magnets stimulated by their interaction with 
the nuclear spin bath. In this model open spins are those 
having the small Zeeman splitting compared to their hy- 
perfine interaction. The sharp kinetic transition is pre- 
dicted for the 2D lattice of magnetic molecules, however 
the abrupt change of the relaxation rate is smeared out 
by its continuous dependence on the total longitudinal 
bias Eq. pip , which may obscure the experimental ob- 
servation of the transition point. 
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VII. APPENDIX 

Consider the probability P{N) that two spins sepa- 
rated by N sites are connected through percolating sites. 
It is clear that if this probability decreases with N slower 
than the inverse number of paths of length N starting at 
the given site, (z— 1)^^, then the infinite cluster of perco- 
lating sites is formed. This criterion is also applicable to 
the static percolation problem on BL where the connec- 
tion probability is given by Pj^. Accordingly, the infinite 
cluster is formed at Pg ^ in full agreement with 

the exact solution^. 

Each site of the line connecting spins and + 1 
has z — 2 neighbors which do not belong to the line. It 
is convenient for the rest of the discussion to introduce 
a short notation for the probability of finding a non- 
percolating spin (or n-spin) when m = 0,1,2 of them 
are definitely percolating, 2 — to neighbors are definitely 
non-percolating and the remaining group of z — 2 neigh- 
bors can contain both percolating and non-percolating 
spins. 



k=0 



{z^2)\{l^WefW^ 

k\{z - fc - 2)! 



-fc-2 



(13) 



cf. Eqs. ([3]) and (jl]). Indeed, P„i are all we need to 
know to calculate the probability that a given line spin is 
non-percolating when the states of its line neighbors are 



known. Since We is known as a function of Wq, we can 
concentrate on properties of the line spins alone. Note, 
that this decomposition of probabilities Eq. ([T3|) is pos- 
sible due to the Bethe lattice factorization around non- 
percolating spins. In other words different realizations 
for spins belonging to branches of BL separated by n- 
spins are independent of each other. Since the states of 
the end spins (0 and A^ -|- 1) can not change scaling prop- 
erties of the long line connecting them we select them 
to be n-spins to simplify the derivation of P{N). We 
then consider P{N) as complementary to the probability 
of having an arbitrary line decomposition in terms of p- 
and n-spins with at least one n-spin among A^ > 0, i.e. 



1 



{"} 



(14) 



where in each sequence {n} all A^ + 2 spins are divided 
into 2s-|-l segments of alternating n- andp-spins as {71} = 
(ni,pi,n2 . . .ps,ns+i) with ni,pi > 1 [ui is used for n- 
spins and pi for p-spins). The probability of a particular 
sequence is denoted as p„^p[ni,pi, . . . ,ps,ns-|_i]. Clearly, 

"1 -I- pi H VPs+ ns+i = N + 2. 

Due to factorization provided by n-spins (no informa- 
tion can be exchanged between the neighbors of non- 
percolating spins and thus all branches of BL around 
them are statistically independent) we have 



Pn.p 



' 



(for s = or ni = A^ + 2) , 



(15) 



iS-l 



Pn,p = Pe(«l)Pe("s+l )X{pn{nj)X{P{p,) , (for s>0) 
J=2 3 = 1 

Here p„(to) is the probability of having a cluster of m 
n-spins in a row, and, Peirn) is a similar quantity for the 
first and last groups of n-spins. These probabilities can 
be immediately computed by counting how many n-spins 
have to p-spins as their neighbors (recall that the end 
spins (0 and A^ + 1) are assumed to be non-percolating) 



Pn(l) =P2, P„(to > 1) = p2j^o™-2 
Pe(l) =1, Peim > 1) = FiF^-\ 



(16) 



The self-consistent equation for P{N) is solved us- 
ing generating functional P{x) = X]^=i a;^P(A^) (and 
similarly defined Pe{x) = X]w=i ^^/'e(-^) ^^'^ Pn{x) = 
Sn=i Pn(N)) which transforms (fTi]) into 



Fox 



pI{x)P{x) 



1 — X 1 — Fox 



^[p„(z)p(z)]^ (17) 



fc=0 



with Pe{x) = :r + FixP' / {1 — Fqx) and pn{x) = F2X 
F^x^ /{I — Fqx). After elementary algebra we find 



P{x) 



a:3(l-Po) 



(1 - .t)(1 - Fox)pI{x) + a;3(l - Fq)p„{x) 



which can be further simplified to 

xil - Fo) 



(18) 



P(x) 



(1 - x)(l - Fox + 2Fix)) + a;2[p2 + ^^(i _ p^)^ 
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The infinite cluster is formed if P{N) decreases with N 
slower than l/(z-l)^). This means that the threshold is 
determined by the divergence of P(x) at a; = (2— 1) which 
is only possible if the denominator in Eq. (|18p is zero. 
Thus the formation of the infinite cluster is determined 
by the equation 

(z-l)(l + Fo -2i^i) = 
1 + - l)'[i^2(l - Fo) + + Fo- 2Fi] , (19) 
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